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Introduction:  Prostate  cancer  (PCa)  is  the  most  common  cancer  in  men  in  several 
countries,  with  the  American  Cancer  Society  (ACS)  estimating  241,740  new  cases  of  PCa 
to  be  diagnosed  and  deaths  of  28,170  men  of  PCa  (1).  Due  to  its  prevalence  in  the  male 
population  as  well  as  its  unpredictable  clinical  course,  early  detection  and  diagnosis  have 
become  a  priority  for  many  health  care  professionals.  Another  method  for  staging 
prostate  cancer  is  through  imaging  techniques  including  ultrasound,  computed 
tomography  (CT),  and  magnetic  resonance  imaging  (MRI)  with  or  without  the  help  of 
dynamic  contrast  enhancement  modeling  (DCE-MRI),  diffusion  weighted  imaging 
(DWI),  and  magnetic  resonance  spectroscopy  (MRS)  (2-5).  MRS  is  a  powerful  tool  for 
exploring  the  cellular  chemistry  of  human  tissues  (3, 5, 6-1 1).  There  is  a  growing  body  of 
evidence  that  1H  MRS  may  contribute  to  the  clinical  evaluation  of  prostate  cancer  and 
also  for  evaluating  the  metabolic  alterations  due  to  therapy.  There  have  been  no  reports 
on  combining  two  spectral  dimensions  with  two-dimensional  (2D)  or  three  dimensional 
(3D)  spatial  encoding  applicable  to  prostate  cancer.  Acceleration  of  magnetic  resonance 
spectroscopic  imaging  (MRSI)  has  been  demonstrated  using  echo-planar  imaging 
techniques  (12-13).  Recently,  Schulte  et  al.  have  successfully  developed  an  algorithm 
called  prior-knowledge  fitting  (ProFit)  to  quantify  metabolite  concentrations  using  the 
JPRESS  spectra  recorded  using  a  Philips  3T  MRI  scanner  (14).  It  was  demonstrated  that 
metabolite  quantitation  of  JPRESS  spectra  with  ProFit  was  accurate,  robust  and  yielding 
generally  consistent  results,  both  in  vivo  and  in  vitro.  Their  results  suggest  that  the 
number  of  quantifiable  prostate  metabolites  can  be  increased  from  3-4  with  ID 
PRESS/LC-Model  to  more  than  10  with  JPRESS/ProFit  (15-16). 

Body: 

i)  Proposed  Task  1  (Months  1-6):  To  implement  a  multi-voxel  based  extension  of  the 
JPRESS  sequence,  in  which  two  spectral  encodings  will  be  combined  with  two  spatial 
encodings  using  the  new  Siemens  VB17a  platform.  This  four-dimensioned  (4D)  data 
acquisition  scheme  will  be  accomplished  utilizing  the  EPI  approach  that  is  commonly 
used  for  spatied  encoding  in  MRI. 

Completed  and  Reported  in  the  1st  year  Annual  Report. 

ii)  Proposed  Task  2:  To  eveduate  the  EPI-based  JPRESS  using  a  prostate  phantom 
containing  several  metabolites  which  have  been  reported  in  prostate  tissues,  and  to 
optimize  the  EP-JJRESI  sequence  and  other  acquisition  parameters  using  the  phantom 

(Months  6-12). 

Completed  and  Reported  in  the  1st  year  Annual  Report. 

iv)  Proposed  Task  3:  To  record  the  4D  EP-JRESI  spectra  in  the  peripheral,  central  and 
transition  zones  of  healthy  prostates.  (Months  6-18). 

Completed  partially  and  Reported  in  the  1st  year  Annual  Report.  More  work  is  in 
progress. 

Proposed  Task  4:  To  develop,  evaluate  and  optimize  the  prior -knowledge  basis  set 
spectra  using  the  GAMMA- simulation  and  prostate  phantom  solutions  as  prior 
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knowledge  for  the  multi-voxel  based  JPRESS  spectra  recorded  using  the  3T  MRl  scanner 

( Months  6-12). 

Accomplished  during  September  29,  2011-October  28  2012: 

Proposed  Task  5:  To  record  the  multi-voxel-based  2D  spectra  in  the  peripheral,  central 
and  transition  zones  of  patients  with  BPH  and  malignant  prostate  cancer.  The  prostate 
metabolite  concentrations  calculated  using  the  ProFit  algorithm  prepared  for  the  multi¬ 
voxel  data  will  be  compared  with  LC-Model  processed  ID  spectral  based  MRSI  data 

(Months  18-36). 

Accomplished  during  October  29,  2012-October  28  2013:  We  present  our  recent 
results  from  2  malignant  prostate  cancer  patients.  The  following  parameters  were  used  to 
acquire  the  4D  EP-JRESI  data:  TR/TE=1.5s/30ms,  2  averages,  512  t2,  oversampled  32kx, 

64  increments  along  the 
indirect  spectral  (ti)  and 
16  spatial  ky  dimensions; 
the  endorectal  "receive" 
coil.  First,  the  4D  NUS 
EP-JRESI  data  was 
acquired  in  a  61  year  old 
PCa  patient  having  PSA  of 
9.1  and  two  malignant 
lesions  (GS3+4  in  the 
right  base  and  GS4+3  in 
the  right  mid  regions)  and 
extracted  spectra  from  2 
different  locations  (lml) 
Figure  1.  i)  a  coronal  T2W  slice  showing  the  multi-voxel  EP-JRESI  grids;  extracted 
2D  JPRESS  spectra  of  malignant  voxels  (lml)  from  the  right  mid  (ii)  and  right  base 
(iii)  portions 

after  the  CS  reconstruction  are  shown  in  Figure  1.  The  top  insert  (i)  shows  the  multi¬ 
voxel  grid  overlaid  on  an  T2W  coronal  slice;  the  extracted  two  spectra  from  the  malignant 
masses  in  the  right  mid  (GS4+3)  (ii)  and  right  base  (iii)  regions. 

Shown  in  the  extracted  tumor  spectra  of  this  61  year  old  patient  with  two  different 
lesions  (GS3+4  and  GS  4+3)  are  the  following:  significantly  increased  Ch  and 
significantly  decreased  Cit,  Spin,  ml,  Tau  and  si  demonstrating  more  metabolite  changes 
than  the  commonly  reported  few  metabolite  changes  using  the  conventional  3D  MRSI 
technique  with  average  weighted  encoding.  It  is  also  seen  that  the  lesion  with  GS4+3 
shows  more  significant  changes  of  above  mentioned  metabolites,  and  also,  elevated  Glx 
(Glu+Gln)  in  agreement  with  a  recent  HR-MAS  analysis  of  radically  resected  prostate 
specimens  by  Gribbestad  and  co-workers  (17).  Presence  of  an  additional  peak  in  GS4+3 
at  2.8ppm  was  indicative  of  the  elevated  to 6  fatty  acids  as  reported  recently  (18) 

Second,  shown  in  Fig. 2  is  the  4D  NUS  EP-JRESI  data  acquired  in  a  72  year  old  PCa 
patient  having  PSA  of  5.6  and  GS  4+3  and  extracted  spectra  from  a  cancer  location 
location  (lml). 
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Figure  2.  i)  An  axial  ToW  MRI  slice  showing  the  multi-voxel  EP-JRES1  grids;  ii)  multi¬ 
voxel  display  of  Ch  and  Cr  peaks  and  an  extracted  2D  JPRESS  spectrum  from  one  of 
malignant  voxels  (iii). 

Key  Research  Accomplishments 

•  Continued  evaluation  of  the  4D  EP-JRESI  sequence  in  6  malignant  prostate 
cancer  patients  and  one  healthy  male  subject. 

•  We  have  made  further  progress  on  compressed  sensing  reconstruction  of  the 
non-unifonnly  undersampled  4D  EP-JRESI  sequence  using  different 
reconstruction  methods  such  as  maximum  entropy,  total  variation  (TV),  etc  .  This 
will  facilitate  shortening  the  endorectal  spectral  acquisition  and  reducing  the 
patient  inconvenience  during  the  scan. 

•  Using  the  preliminary  results  obtained  using  this  IDEA  grant,  we  had  reported  our 
failed  attempt  on  an  R01  application  entitled  "EchoPlanar  J-resolved  Prostate 
Cancer  Metabolite  Imaging  Using  Compressed  Sensing"  that  was  submitted  to  the 
National  Institute  of  health  (NIH)  in  July20 12/October  2011  (32nd  percentile). 
Hence,  a  new  NIH  R01  grant  application  entitled  “Fast  J-resolved  Prostate  MR 
Spectroscopic  Imaging  and  Non-linear  Reconstruction”  was  submitted  on  October 
5,  2013  using  the  preliminary  results  from  this  IDEA  grant. 
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Reportable  Outcomes: 

A.  Peer-reviewed  Publications: 

1.  Nagarajan  R,  Margolis  D,  Raman  S,  et  al.  MR  Spectroscopic  Imaging  of  Peripheral 
zone  in  Prostate  Cancer  using  a  3T  MRI  scanner:  Endorectal  versus  External  Phased 
array  Coils.  Magnetic  Resonance  Insights  2013:  6:  51-58  [PMID:  N/A]. 

2.  Thomas  MA,  Nagarajan  R,  Huda  A,  et  al.  Multidimensional  MR  Spectroscopic 
imaging  of  Prostate  Cancer  In  vivo.  NMR  in  Biomed  2013  July  3  l(Epub  ahead  of  print). 
[PMID:  23904127]. 


B.  Presentations:  During  the  2nd  year,  the  first  abstract  entitled  “Accelerated  Four- 
Dimensional  Echo-Planar  J-Resolved  Spectroscopic  Imaging  of  Human  Prostate: 
Prospective  Non-Uniform  Undersampling  and  Maximum  Entropy  Reconstruction.”  was 
submitted  and  presented  at  the  54th  Experimental  Nuclear  Magnetic  Resonance  (ENC), 
Asilomar,  CA,  Abstract  #214,  (April  14-19,  2013).  The  2nd  abstract  entitled 
“Nonuniformly  Under-Sampled  (NUS)  Echo  Planar  J-Resolved  Spectroscopic  Imaging 
(EP-JRESI)  of  Prostate  Cancer  Patients  and  Compressed  Sensing  Reconstruction” 
summarizing  the  implementation  of  the  4D  EP-JRESI  sequence  and  evaluation  of  it  in  the 
prostate  cancer  patients  was  presented  at  the  21st  International  Society  of  Magnetic 
Resonance  in  Medicine  (ISMRM),  Salt  Lake,  Utah,  April  20-26,  2013,  #3959. 

C.  Books:  None  on  Prostate  Cancer  Research  based. 

Conclusions:  After  successfully  implementing  the  4D  EP-JRESI  scanning  protocol  on 
the  3T  MRI  scanner  and  testing  it  in  healthy  controls  and  patients  during  the  1st  year,  7 
more  subjects  have  been  scanned  in  this  year.  We  will  continue  to  recruit  20  malignant 
and  5  benign  prostate  cancer  patients,  and  5  healthy  male  subjects  during  the  next  year. 
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Abstract:  Magnetic  resonance  spectroscopic  imaging  (MRSI)  detects  alterations  in  major  prostate  metabolites,  such  as  citrate  (Cit), 
creatine  (Cr),  and  choline  (Ch).  We  evaluated  the  sensitivity  and  accuracy  of  three-dimensional  MRSI  of  prostate  using  an  endorectal 
compared  to  an  external  phased  array  “receive”  coil  on  a  3T  MRI  scanner.  Eighteen  patients  with  prostate  cancer  (PCa)  who  under¬ 
went  endorectal  MR  imaging  and  proton  (1H)  MRSI  were  included  in  this  study.  Immediately  after  the  endorectal  MRSI  scan,  the  PCa 
patients  were  scanned  with  the  external  phased  array  coil.  The  endorectal  coil-detected  metabolite  ratio  [(Ch+Cr)/Cit]  was  significantly 
higher  in  cancer  locations  (1.667  +  0.663)  compared  to  non-cancer  locations  (0.978  ±  0.420)  (P  <  0.001).  Similarly,  for  the  external 
phased  array,  the  ratio  was  significantly  higher  in  cancer  locations  (1.070  ±  0.525)  compared  to  non-cancer  locations  (0.521  ±  0.310) 
(P  <  0.001).  The  sensitivity  and  accuracy  of  cancer  detection  were  81%  and  78%  using  the  endorectal  ‘receive’  coil,  and  69%  and  75%, 
respectively  using  the  external  phased  array  ‘receive’  coil. 

Keywords:  prostate  cancer,  MRSI,  metabolites,  external  coil,  sensitivity,  citrate,  choline,  creatine 
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Introduction 

Prostate  cancer  (PCa)  is  the  most  common  non- 
cutaneous  cancer  and  second  leading  cause  of  cancer 
death  in  men.  In  2012,  approximately  241,740  new 
cases  and  28,170  PCa-related  deaths  occurred  in  the 
United  States. 1  Magnetic  resonance  imaging  (MRI)  has 
been  used  to  evaluate  prostate  anatomy  and  prostate 
pathologies  for  several  years.  MRI,  with  its  excellent 
soft-tissue  differentiation,  provides  high-resolution 
images  of  the  prostate  and  surrounding  structures. 
Magnetic  resonance  spectroscopy  (MRS)  is  a  power¬ 
ful  tool  for  exploring  the  cellular  chemistry  of  human 
tissues.2-6  There  is  a  growing  body  of  evidence  that 
proton  MRS  can  be  used  for  clinical  evaluation  of  PCa 
and  metabolic  alterations  before  and  after  therapy. 
Magnetic  resonance  spectroscopic  imaging  (MRSI) 
can  be  used  to  measure  metabolite  levels  in  the  tissue, 
particularly  choline  (Ch),  citrate  (Cit),  creatine  (Cr), 
and  various  polyamines  (spennine,  spermidine,  and 
putrecine).  PCa  typically  shows  an  increased  concen¬ 
tration  of  Ch  and  reduction  of  Cit  and  polyamines. 

MRSI  of  the  prostate  is  typically  perfonned  using  a 
combination  of  point-resolved  spectroscopy  (PRESS)7 
and  three-dimensional  localized  MRSI8  rather  than 
the  traditional  single-voxel  or  slice-based  two-dimen¬ 
sional  MRSI  technique  was  used  in  brain  spectro¬ 
scopic  imaging.  3D  MRSI  requires  phase  encoding 
along  the  three  spatial  dimensions,  conventionally 
known  as  frequency,  phase,  and  slice.  Acquisition 
time  and  coverage  of  the  prostate  are  the  primary  con¬ 
siderations  in  choosing  matrix  dimensions.  Acquiring 
3D  MRSI  data  with  higher  spatial  resolution  requires 
a  long  total  acquisition  time.  Conventional  prostate 
MRSI  studies  involving  average  weighted  encoding 
use  a  long  echo  time  (TE)  with  a  short  repetition  time 
(TR),  allowing  observation  of  a  reduced  number  of 
metabolites.9  Long  TEs  are  used  for  3D  MRSI  due 
to  the  addition  of  MEGA10  radio-frequency  pulses 
for  both  water  and  lipid  suppression.  The  interpreta¬ 
tion  system  most  used  to  discriminate  between  cancer 
and  normal  prostatic  tissue  in  the  peripheral  zone  was 
described  by  Kurhanewicz  et  al.9  They  calculated  the 
peak  area  ratios  of  Ch  and  Cr  to  Cit  [(Ch+Cr)/Cit]  for 
each  voxel.  Inclusion  of  Cr  in  this  ratio  is  mandatory 
because  of  the  Cr  peak  is  very  close  to  the  Ch  peak  in 
the  spectrum.  However,  Cr  appears  be  maintained  at 
a  relatively  constant  level  in  both  healthy  and  tumor 
prostatic  tissues. 


It  was  unknown  whether  the  quality  of  endorectal 
MRSI  scanning  could  be  improved  by  using  a  stron¬ 
ger  (3T)  scanner  and/or  using  perfluorocarbon  (PFC) 
in  the  endorectal  coil  (ERC)  instead  of  air.  Rather 
than  filling  air,  PFC  has  been  shown  to  improve 
MRSI  image  quality  and  the  magnetic  susceptibil¬ 
ity  closely  matching  that  of  prostate.11  The  ERC  is 
contraindicated  in  a  number  of  patients  (eg,  after 
abdominal  perineal  resection  for  rectal  cancers  or 
after  radiation  therapy  for  the  pelvis).12  In  some  cases, 
it  is  preferable  to  not  use  an  ERC  to  avoid  structural 
deformation  of  the  prostate  peripheral  zone,  which  is 
often  compressed  by  an  ERC.  The  ERC  also  causes 
signal  hyperintensity  near  the  rectum  and  the  neigh¬ 
boring  peripheral  zone.  Signal  hyperintensity  and 
tissue  deformation  can  make  complicate  diagnostic 
interpretation.  Previous  studies  have  reported  the  use 
of  external  phased  array  coils  for  prostate  MR  spec¬ 
troscopy  in  3  T. 13-20  The  major  goal  of  the  study  was  to 
evaluate  3T  MRSI  for  the  peripheral  zone  in  prostate 
cancer  patients  using  an  endorectal  coil  and  to  com¬ 
pare  the  performance  of  3T  MRSI  and  an  external 
phased  array  ‘receive’  coil. 

Materials  and  Methods 

Eighteen  patients  ranging  in  age  from  56-72  years 
(mean,  63.1  years)  with  PCa  who  underwent  endorec¬ 
tal  MR  imaging  and  proton  MR  spectroscopic  imag¬ 
ing  were  included  in  this  study.  The  study  protocol 
was  approved  by  the  Institutional  Review  Board,  and 
informed  consent  was  obtained  from  each  patient. 
Gleason  scores  for  the  tumor  on  prostate  biopsy 
ranged  from  3  to  8,  while  prostate-specific  antigen 
varied  from  2.8  to  20.6  ng/mL  (mean  of  6.84  ng/mL). 
A  Siemens  3T  MRI  Scanner  with  high-performance 
gradients  (Trio-Tim,  Siemens  Medical  Solutions, 
Erlangen,  Germany)  was  used  in  this  investigation. 
A  quadrature  body  ‘transmit’  coil  was  used  to  trans¬ 
mit  radio-frequency  pulses.  An  endorectal  inflatable 
‘receive’  coil  (Medrad  Corporation,  Indianola,  PA, 
USA)  was  then  inserted  into  the  rectum  and  inflated 
with  50  cc  PFC.  The  coil  was  positioned  horizon¬ 
tally  at  approximately  the  10:00  and  2:00  positions. 
After  endorectal  scanning,  patients  were  scanned  with 
the  external  phased  array  ‘receive’  coil  for  the  com¬ 
parison  study,  ff-wcightcd  images  in  the  transverse, 
sagittal,  and  coronal  planes  were  acquired  by  using  a 
turbo  spin-echo  sequence.  The  MRI  protocol  included 
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T,-wcightcd  images  acquired  using  a  fast  spin-echo 
sequence  with:  repetition  time  (TR)  =  3,800  ms, 
effective  echo  time  (TE)  =  101  ms,  slice  thick¬ 
ness  =  3  mm,  field  of  view  =140  mm,  and  matrix 
size  =  256  x  256  mm2.  MRSI  was  performed  in  all 
patients,  which  included  a  3D  water-  and  fat-suppressed 
spectroscopic  acquisition.  3D  MRSI  parameters  of  the 
endorectal  and  external  phased  array  coil  were  as  fol¬ 
lows:  TR  750  ms,  TE  145  ms,  acquisition  bandwidth 
1250  Hz,  6  averages,  and  5 12  spectral  data  points  with 
a  voxel  resolution  of  0.3  mL.  Total  acquisition  time 
was  approximately  12  min.  For  the  external  phased 
array,  the  voxel  resolution  was  0.35  mL.  A  PRESS- 
based  sequence  was  used  to  acquire  proton  MR  spec¬ 
tra  from  a  volume  of  interest  (VOI)  of  approximately 
55  x  45  x  45  mm3.  Outer  volume  suppression  of  water 
and  lipid  was  achieved  using  eight  3 -cm  thick  satura¬ 
tion  pulses  around  the  VOI. 

A  spectroscopist  examined  the  MRSI  data  set 
and  reported  the  location  and  number  of  suspicious 
voxels  to  the  radiologist,  to  identify  metabolite 
ratios  predictive  of  cancer.  For  3D  MRSI  post-pro¬ 
cessing,  each  spectrum  was  Fourier  transformed, 
frequency-,  phase-,  and  baseline-corrected,  and  the 
peaks  of  Cit,  Ch,  and  Cr  were  subsequently  fitted. 
A  Hamming  filter  was  used  for  the  MRSI  spatial 
dimensions  of  the  data.  For  tumor  localization,  the 
prostate  was  split  along  the  midline  and  further 
divided  into  the  apex,  the  middle,  and  the  base  of 
the  gland.  According  to  histopathologic  findings, 
the  voxels  and  corresponding  spectra  were  assigned 
as  normal  tissue  or  tumor  tissue.  For  selected  vox¬ 
els  in  the  peripheral  zone,  the  area  under  the  curve 
of  the  metabolite  resonances  was  determined  and 
the  signal  intensity  ratio  for  (Ch+Cr)/Cit  was  calcu¬ 
lated  using  commercially  available  software  from 
Siemens.  Peak  areas  for  Ch,  Cr,  Spin,  and  Cit  were 
calculated  using  numeric  integration.  Metabolic 
maps  of  (Ch+Cr)/Cit  were  generated  since  Spin 
cannot  be  separated  from  Ch  and  Cr  peaks.  Voxels 
were  considered  suitable  if  they  consisted  of  at  least 
75%  peripheral  zone  tissue  and  did  not  include  peri¬ 
urethral  tissue.  If  the  Cit  peak  was  lower  than  the 
Ch  peak  or  was  undetectable,  the  voxel  was  deter¬ 
mined  to  be  malignant.  If  the  Cit  peak  was  higher 
than  the  Ch  peak,  the  voxel  was  considered  noncan- 
cerous  for  (Ch+Cr)/Cit  values  smaller  than  0.50  and 
malignant  for  (Ch+Cr)/Cit  values  greater  than  0.50. 


For  the  external  phased  array  coil,  the  noncancer- 
ous  ratio  was  less  than  0.28  and  the  malignant  ratio 
was  greater  than  0.28. 

The  paired  /-test  was  used  to  determine  whether 
the  ratios  of  (Ch+Cr)/Cit  in  cancer  were  different 
from  those  in  non-cancer  using  an  endorectal  coil  and 
external  coil.  AR-value  less  than  0.05  was  considered 
statistically  significant. 

Receiver  operating  characteristic  (ROC)  curve 
analyses  based  on  logistic  regression  models  were 
perfonned  to  identify  the  optimal  cutoff  value  for  pre¬ 
dicting  metabolite  ratios  using  the  endorectal  coil  and 
external  phased  array  coil.  The  area  under  the  curve 
(AUC),  interpreted  as  the  average  value  of  sensitiv¬ 
ity  for  all  possible  values  of  specificity,  was  used  in 
the  ROC  analysis.  An  area  of  0.50  implies  that  the 
variable  adds  no  information,  whereas  an  area  of  1 
implies  perfect  accuracy.  Sensitivity,  specificity,  posi¬ 
tive  predictive  value  (PPV),  and  negative  predictive 
value  (NPV)  as  well  as  accuracy  were  reported  for  the 
optimal  thresholds.  P- values  <0.05  were  considered 
statistically  significant.  Statistical  analyses  were  per¬ 
fonned  using  SPSS  software  (Statistical  Package  for 
the  Social  Sciences,  Version  18.0,  SPSS  Inc,  Chicago, 
IL,  USA). 

Results 

Figure  1 A  and  B  show  the  metabolite  ratios  (Ch+Cr)/ 
Cit  for  cancer  and  non-cancer  locations  with 
endorectal  and  external  phased  array  coil  in  18  PCa 
patients.  For  the  endorectal  coil,  the  metabolite  ratio 
(mean  ±  SD)  was  significantly  higher  in  cancer  loca¬ 
tions  (1 .667  ±  0.663)  compared  to  in  non-cancer  loca¬ 
tions  (0.978  ±  0.420)  (P  <  0.001).  Similarly,  for  the 
external  phased  array  coil,  the  ratio  was  significantly 
higher  in  cancer  locations  (1.070  ±  0.525)  than  in 
non-cancer  locations  (0.521  ±  0.310)  ( P  <  0.001). 
Figure  2 A  and  B  show  a  comparison  of  endorectal 
and  external  phased  array  MR  Spectroscopic  Imag¬ 
ing  of  a  67  year-old  PCa  patient.  A  significant  eleva¬ 
tion  in  Ch/Cr  and  decreased  Cit  were  observed  in  the 
right  side  peripheral  zone  of  the  PCa  patient  scanned 
using  endorectal  coil.  Similar  trends  for  metabolite 
changes  are  shown  in  Figure  2B  with  a  slightly  worse 
signal-to-noise  ratio  (SNR)  because  the  patient  was 
scanned  using  an  external  coil.  Table  1  shows  the  full 
width  at  half  maximum  (FWHM)  of  Cit  and  (Cr+Ch) 
in  cancer  and  noncancerous  locations  using  the 
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Endorectal  patients 


External  coil  patients 


Figure  1.  Comparison  of  metabolite  ratios  [(Ch+Cr)/Cit]  in  18  PCa  patients  scanned  with  the  endorectal  and  external  body  array  ‘receive’  coils. 
Note:  Significant  (P  <  0.05)  elevation  of  Ch/Cr  observed  in  cancer  locations  compared  to  non-cancer  locations  in  both  coils. 


endorectal  and  external  phased  array  coil.  The  mean 
FWHM  of  Cit  in  cancer  and  non-cancer  locations  for 
patients  scanned  using  the  endorectal  coil  were  11.7 
and  11.5  Hz  respectively,  whereas  using  the  external 
coil,  the  mean  FWHM  values  were  17.7  and  19.5  Hz, 
which  showed  a  higher  coefficient  of  variance  (CV). 
Similarly,  the  mean  FWHM  of  Ch+Cr  in  cancer  and 
non-cancer  locations  of  patients  scanned  using  the 
endorectal  coil  were  17.7  and  19.5  Hz,  where  in  the 
external  coil,  the  mean  FWHM  values  were  25.7  and 
28  Hz  with  a  higher  CV.  The  mean  global  water  line 
width  of  25  Hz  observed  using  the  endorectal  coil  and 


3 1  Hz  was  observed  using  the  external  phased  array 
coil  with  higher  CV. 

Figure  3  shows  that  the  MRSI  metabolite  ratios 
in  peripheral  zone  cancer  locations  derived  from 
endorectal  coil-  and  external  body  coil  were  positively 
and  linearly  correlated  (R2  =  0.571).  Figure  4  shows 
the  comparison  of  receiver  operating  characteristics 
curves  for  the  endorectal  coil-  and  external  body  coil- 
derived  MRSI  ratio.  ROC  curve  analyses  for  differ¬ 
entiating  endorectal  coil  suggested  an  optimal  cutoff 
value  of  1 .35.  This  implies  that  the  proportions  of  cor¬ 
rectly  identified  endorectal  coil  sensitivity,  specificity, 


54 


Magnetic  Resonance  Insights  2013:6 


MRSI  of  prostate  cancer  using  endorectal  and  external  coils 


Ch/Spm/Cr  Ch/Spm/Cr 


Figure  2.  Comparison  of  (A)  endorectal  and  (B)  external  body  array  MR  Spectroscopic  Imaging  of  a  67-year  old  PCa  patient. 

Note:  Extracted  single  voxel  spectra  show  a  significant  (P  <  0.05)  elevation  of  Ch/Cr  and  declined  Cit  in  the  right  side  peripheral  zone  of  the  PCa  patient 
in  both  coils. 


PPV,  and  NPV  were  8 1 .3%,  75.0%,  76.5%,  and  80.0%, 
respectively.  The  AUC  was  86.9%,  with  an  accuracy 
of  78.1%.  The  cutoff  value  for  the  external  body  coil 
was  0.77  for  sensitivity,  specificity,  and  accuracy 
values  of  68.8%,  81.3%,  and  75.0%,  respectively. 
The  AUC  of  the  MRSI  ratio  slightly  increased  for 
endorectal  coil  (86.9%)  compared  to  external  phased 
array  coil  (85.9%).  In  the  evaluation  of  the  endorec¬ 
tal  coil  versus  external  phased  array  coil,  the  MRSI 
ratios  showed  an  accuracy  of  78.1%,  indicating  good 
discrimination  compared  to  external  phased  array  of 
75.0%.  Detailed  results  of  the  ROC  curve  analyses 


are  shown  in  Table  2,  which  shows  the  sensitivity, 
specificity,  PPV,  NPV,  AUC,  and  accuracy  of  classi¬ 
fying  endorectal  coil  and  external  body  coil. 

Discussion 

Numerous  studies  have  suggested  that  MRSI  using 
an  ERC  is  the  most  promising  technique  for  detect¬ 
ing  and  staging  of  PCa.  At  more  commonly  available 
clinical  field  strengths  of  1.5T,  an  ERC  is  necessary 
for  obtaining  a  sufficiently  high  SNR  with  subsequent 
spatial  resolution,  allowing  reliable  cancer  delineation 
in  a  clinically  reasonable  time  frame.  However,  the 


Table  1.  FWHM  of  water  and  metabolite  resonances  using  endorectal  and  external  body  array  “receive”  coils. 


Peak 

Endorectal  FWHM 
Mean  ±  SD 

CV  (%) 

External  body  array  FWHM 
Mean  ±  SD 

CV  (%) 

Cancer  Cit 

11.73  ±2.75 

23.41 

15.45  ±  6.17 

39.96 

Cancer  (Cr+Ch) 

17.71  ±6.78 

27.57 

25.75  ±8.47 

32.88 

Non-cancer  Cit 

11.53  ±  1.73 

15.02 

1 8.83  ±  10.16 

53.94 

Non-cancer  (Cr+Ch) 

19.59  ±6.83 

27.77 

28.06  ±  10.80 

38.49 

Water 

25.00  ±5.16 

20.64 

31.00  ±9.62 

31.05 

Note:  The  FWHM  of  the  water  resonance  represents  global  line  width  whereas  the  FWHM  values  of  metabolites  were  derived  from  spectra  extracted  from 
selected  locations. 
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Figure  3.  Correlation  between  the  maximum  ratios  [(Ch+Cr)/Cit]  of  can¬ 
cer  identified  in  the  peripheral  zone  with  external  phased  array  coil  vs. 
endorectal  coil  MRS. 

Note:  A  positive  correlation  was  observed. 


use  of  an  ERC  is  more  time-consuming,  leading  to 
higher  costs  and  greater  patient  discomfort.  Thus,  we 
compared  the  performance  of  the  endorectal  coil  with 
the  external  phased  array  coil  at  3T. 

Compared  with  external  phased  array  coil  at  3T, 
endorectal  MR  spectroscopy  at  3T  significantly 
improves  spectral  line  width  and  coefficient  of  vari¬ 
ance  of  metabolite  ratio.  The  FWHM  is  a  good  indi¬ 
cator  of  spectral  homogeneity.  The  FWHM  value 
typically  increases  as  VOI  increases.  Excellent  quality 
of  prostate  spectra  can  be  obtained  through  optimiza¬ 
tion  of  B0  homogeneity  for  the  PRESS-selected  VOI. 
This  procedure  typically  involves  the  combined  use 
of  a  standard  automatic  shim  provided  by  the  man¬ 
ufacturer  and,  if  necessary,  manual  adjustments  of 
the  linear  x,  y,  and  z  gradients.  Small  improvements 
in  shim  can  make  a  significant  difference  in  spectra 
quality.  Particularly,  good  B0  homogeneity  is  essential 


for  sufficient  water  and  lipid  suppression.  Water  and 
lipid  suppression  is  achieved  through  generation  of 
frequency-selective  MEGA  pulses.10  During  MRSI 
acquisition,  large  lipid  resonance  peaks  obscures  the 
metabolite  peak.  Outer  volume  saturation  pulses  were 
used  to  eliminate  signals  from  adjacent  tissues,  par¬ 
ticularly  periprostatic  lipids  and  rectal  wall  tissue. 
The  coefficient  of  variance  was  higher  in  the  external 
phased  array  than  for  the  endorectal  coil  due  to  the 
larger  size  of  the  coil  and  increased  distance  from  the 
prostate  proportional  to  abdominal  circumference, 
leading  to  a  decrease  in  the  SNR  ratio. 

The  overall  accuracy  of  the  external  phased  array 
coil  MRSI  has  been  shown  to  be  slightly  inferior  to 
that  using  the  endorectal  coil  in  our  pilot  study  involv¬ 
ing  a  limited  number  of  patients.  In  the  localization  of 
PCa  with  MRSI  at  3T,  the  use  of  the  endorectal  coil 
showed  a  significantly  (P  <  0.05)  higher  sensitivity 
(81.3%)  than  the  external  phased  array  (68.8%).  In 
contrast,  the  specificities  of  cancer  detection  were 
81%  and  75%  using  the  external  phased  array  and  the 
endorectal  coils,  respectively.  Additionally,  there  was 
a  slight  increase  in  the  AUC  (86.9%)  for  endorec¬ 
tal  compared  to  the  external  phased  array  AUC 
(85.9%). 

In  our  investigations  using  an  external  phased 
array  coil,  the  effective  voxel  size  was  slightly 
larger  to  compensate  for  the  poor  SNR  compared 
to  examinations  of  the  prostate  using  an  endorectal 
coil.  Overlap  exists  in  the  (Ch+Cr)/Cit  ratios  between 
cancer  and  noncancerous  tissues.  For  both  ERC  and 
external  phased  array  coils,  (Ch+Cr)/Cit  ratios  were 
significantly  higher  in  cancer  locations  than  non¬ 
cancer  locations  in  this  study,  in  accordance  with  the 


0.8 

> 

!  o.6 

'35 

c  0  4 

<D 

W 

0.2 

0 


Figure  4.  ROC  curves  of  MRSI  ratios  of  prostate  cancer  patients  (A)  with  the  endorectal  and  (B)  the  external  phased  array  ‘receive’  coils. 
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Table  2.  Measures  of  sensitivity,  specificity,  PPV,  NPV, 
and  accuracy  of  endorectal  and  external  body  coil  MRSI 
ratios  in  discrimination  of  cancer  from  non-cancer  with 
ROC  curve  analysis. 


Parameter 

Endorectal  coil 

External  phased  array 

CV 

1.35 

0.77 

Sensitivity,  % 

81.3 

68.8 

Specificity,  % 

75.0 

81.3 

PPV% 

76.5 

78.6 

NPV% 

80.0 

72.2 

AUC,  % 

86.9 

85.9 

Accuracy,  % 

78.1 

75.0 

Abbreviations:  CV,  cutoff  value;  PPV,  positive  predictive  value;  NPV, 
negative  predictive  value;  AUC,  area  under  the  curve. 


results  of  previous  studies.21-30  Additionally,  the  trend 
of  accuracy  and  sensitivity  of  the  ERC  in  this  study 
agreed  with  that  of  a  previous  study.20  The  metabo¬ 
lite  mean  ratio  of  the  external  phased  array  coil  was 
lower  in  cancer  and  non-cancer  locations  due  to  the 
low  SNR.  However,  the  ratio  of  (Ch+Cr)/Cit  using 
both  coils  was  generally  positively  and  linearly  cor¬ 
related.  Although  AUC  and  sensitivity  were  elevated 
for  ERC,  the  results  imply  that  both  ERC  and  exter¬ 
nal  phased-array  coil  MRS  are  feasible  for  detecting 
the  peripheral  zone  cancer  with  an  optimal  threshold 
ratio  for  (Ch+Cr)/Cit. 

Although  the  sensitivity  and  accuracy  increased 
in  ERC,  there  are  some  disadvantages  of  using  this 
method.  Disadvantages  of  using  an  ERC  include  its 
invasiveness,  image  distortions,  and  a  limited  field 
of  view.  The  accuracy  of  endorectal  coil  MRI  is  also 
frequently  affected  by  image  degradation.  The  three 
most  common  causes  of  this  are  near-field  endorectal 
surface  coil  profile  (because  the  coil  is  close  to  the 
gland),  phase  encoding  artifact  (from  motion  of  feces 
or  the  rectum  itself),  and  incorrect  placement  of  the 
coil  (where  the  coil  is  not  parallel  to  the  transverse 
plane  of  the  gland).31 

Another  limitation  is  the  small  number  of  patients 
in  this  study.  Additionally,  there  was  no  correlation 
between  Gleason  scores  and  MRSI  ratios.  Another 
limitation  of  this  study  was  the  inclusion  of  only  the 
peripheral  zone  of  the  prostate,  in  which  75%  cancer 
arises. 

Conclusion 

These  preliminary  findings  confirmed  that  the  use 
of  ERC  significantly  improves  spectral  line  width 


and  the  coefficient  of  variance  of  metabolite  ratio 
compared  with  the  external  phased  array  coil. 
However,  based  on  overall  performance,  use  of  the 
external  phased  array  coil  may  be  recommended  for 
patients  with  rectal  diseases  or  patients  who  could 
not  tolerate  the  discomfort  of  endorectal  surface 
coil  insertion. 
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Prostate  cancer  (PCa)  is  the  second  most  common  type  of  cancer  among  men  in  the  United  States.  A  major  limitation 
in  the  management  of  PCa  is  an  inability  to  distinguish,  early  on,  cancers  that  will  progress  and  become  life 
threatening.  One-dimensional  (ID)  proton  (1H)  MRS  of  the  prostate  provides  metabolic  information  such  as  levels 
of  choline  (Ch),  creatine  (Cr),  citrate  (Cit),  and  spermine  (Spm)  that  can  be  used  to  detect  and  diagnose  PCa.  Ex  vivo 
high-resolution  magic  angle  spinning  (HR-MAS)  of  PCa  specimens  has  revealed  detection  of  more  metabolites  such 
as  myo-inositol  (ml),  glutamate  (Glu),  and  glutamine  (Gin).  Due  to  the  J-modulation  and  signal  overlap,  it  is  difficult 
to  quantitate  Spm  and  other  resonances  in  the  prostate  clearly  by  single-  and  multivoxel-based  1 D  MR  spectroscopy. 
This  limitation  can  be  minimized  by  adding  at  least  one  more  spectral  dimension  by  which  resonances  can  be  spread 
apart,  thereby  increasing  the  spectral  dispersion.  However,  recording  of  multivoxel-based  two-dimensional  (2D) 
MRS  such  as  J-resolved  spectroscopy  (JPRESS)  and  correlated  spectroscopy  (L-COSY)  combined  with  2D  or 
three-dimensional  (3D)  magnetic  resonance  spectroscopic  imaging  (MRSI)  using  conventional  phase-encoding  can 
be  prohibitively  long  to  be  included  in  a  clinical  protocol.  To  reduce  the  long  acquisition  time  required  for  spatial 
encoding,  the  echo-planar  spectroscopic  imaging  (EPSI)  technique  has  been  combined  with  correlated  spectroscopy 
to  give  four-dimensional  (4D)  echo-planar  correlated  spectroscopic  imaging  (EP-COSI)  as  well  as  J-resolved 
spectroscopic  imaging  (EP-JRESI)  and  the  multi-echo  (ME)  variants.  Further  acceleration  can  be  achieved  using 
non-uniform  undersampling  (NUS)  and  reconstruction  using  compressed  sensing  (CS).  Earlier  versions  of  2D  MRS, 
theory  of  2D  MRS,  spectral  apodization  filters,  newer  developments  and  the  potential  role  of  multidimensional 
MRS  in  PCa  detection  and  management  will  be  reviewed  here.  Copyright  ©  2013  John  Wiley  &  Sons,  Ltd. 
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spectroscopic  imaging 


INTRODUCTION 

Prostate  cancer  (PCa)  is  the  second  leading  cause  of  cancer 
deaths  after  lung  cancer  (1).  The  most  common  diagnostic  tools 
used  to  look  for  evidence  of  PCa  include  digital  rectal 
examination  (DRE),  serum  concentration  of  prostate  specific 
antigen  (PSA)  and  transrectal  ultrasonography  (TRUS)  guided 
prostate  biopsy  (2).  The  primary  system  for  prognosis  of  PCa  is 
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based  on  the  size  and  volume  of  cancer  in  a  prostate  biopsy 
specimen  graded  microscopically  according  to  the  glandular 
pattern  of  cancers  (Gleason  score,  GS)  (3).  The  GS  is  the  sum  of 
the  two  most  common  patterns  (grades  1-5)  of  tumor  growth 
found.  In  needle  biopsy,  it  is  now  recommended  that  the  worst 
grade  should  always  be  included  even  if  present  in  less  than 
5%  (4).  Most  PCa  lesions  are  located  in  the  peripheral  zone  (PZ) 
of  the  prostate  and  may  be  detected  by  DRE  when  the  volume 
is  about  0.2  ml  or  larger  (5-7). 

Digital  rectal  examination  has  a  low  overall  sensitivity  (37%) 
and  low  positive  predictive  value  (PPV)  when  lower  PSA  levels 


Abbreviations  used:  JPRESS,  J-resolved  spectroscopy:  MRSI,  magnetic  reso¬ 
nance  spectroscopic  imaging:  CS,  compressed  sensing;  NUS,  non-uniform 
undersampling;  L-COSY,  localized  correlated  spectroscopy;  PCa,  prostate  can¬ 
cer;  BPH,  benign  prostatic  hyperplasia;  PZ,  peripheral  zone;  PSA,  prostate  spe¬ 
cific  antigen;  Ch,  choline;  Cr,  creatine;  Cit,  citrate;  Spm,  spermine;  Glu, 
glutamate;  Gin,  glutamine;  Tau,  taurine ;  GPC,  glycerylphosphocholine ;  GSH, 
glutathione;  TRUS,  transrectal  ultrasound;  GS,  Gleason  score;  CT,  coherence 
transfer;  T2W,  T2  weighted;  DWi,  diffusion-weighted  imaging ;  ADC,  apparent 
diffusion  coefficient;  MVD,  microvessel  density;  DCE,  dynamic  contrast  en¬ 
hanced;  ID,  one  dimensional;  2D,  two  dimensional;  3D,  three  dimensional; 
4D,  four  dimensional;  SNR,  signal  to  noise  ratio;  EPI,  echo-planar  imaging ; 
EPSI,  echo-planar  spectroscopic  imaging;  SV,  single  volume;  EP-COSI,  echo- 
planar  correlated  spectroscopic  imaging;  EP-JRESI,  echo-planar  J-resolved 
spectroscopic  imaging;  ME,  multi-echo;  EFT,  fast  Fourier  transform;  Profit,  prior 
knowledge  fitting;  CRLB,  Cramer-Rao  lower  bound. 
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(0-3  ng  /ml)  are  encountered  (8).  PSA  measurement  has  yielded 
higher  detection  rates  than  DRE  (7),  but  its  specificity  is  low  (36%) 
owing  to  false-positive  PSA  elevation  under  benign  circumstances, 
such  as  inflammation  or  benign  prostatic  hyperplasia  (BPH)  (9). 
These  inaccurate  tools  often  lead  to  incorrect  diagnoses,  inaccurate 
risk  assessments,  patient  anxiety,  and  less  optimal  therapy  choices  in 
management  of  the  disease.  Hence,  there  is  a  need  for  improved 
PCa  diagnosis  with  better-detection,  localization,  and  sampling. 
Prostate  cancer  is  the  only  major  solid  organ  malignancy  that  lacks 
an  imaging-based  diagnosis. 

Diagnostic  imaging  methods  currently  in  use,  such  as  computer¬ 
ized  tomography  (CT)  and  ultrasound  (US),  cannot  adequately 
detect  PCa.  MRI  offers  exquisite  anatomical  details,  but  may  suffer 
from  poor  specificity  in  detecting  and  grading  PCa  (10,11).  T2- 
weighted  (T2W)  MRI  delineates  most  larger  foci  of  PCa  as  a  region 
of  low  signal  intensity  surrounded  by  high  signal  intensity  (longer 
T2)  of  normal  PZ  tissue  (12,13).  Although  the  sensitivity  of  T2W 
images  for  tumor  detection  is  high,  specificity  is  suboptimal  (14). 

Low  specificity  may,  however,  lead  to  overtreatment,  cause 
great  anxiety  to  patients,  and  result  in  many  unnecessary 
biopsies  in  patients  with  no  or  low  grade  cancers.  On  the  other 
hand,  if  the  patient  with  high  grade  cancer  inappropriately 
chooses  "active  surveillance"  as  a  management  option, 
malignant  cells  can  metastasize  to  the  other  parts  of  the  body 
before  the  cancer  becomes  clinically  evident.  Thus  more 
sensitive  and  specific  non-invasive  tests  are  needed  to  better 
differentiate  indolent  and  aggressive  PCa. 

Techniques  such  diffusion-weighted  imaging  (DWI)  and 
perfusion  may  improve  the  performance  of  MR.  DWI  is  an  MRI 
technique  that  visualizes  molecular  diffusion,  that  is,  the 
Brownian  motion  of  water  molecules  in  biologic  tissues  (15),  by 
applying  two  equally  sized  diffusion-sensitizing  gradients,  which 
are  characterized  by  their  b-values.  The  mobility  is  then 
quantified  by  calculating  the  apparent  diffusion  coefficient 
(ADC),  which  depends  mainly  on  the  cellularity,  cell  size, 
extracellular  space,  and  temperature.  When  DWI  is  combined 
with  T2W  MRI,  both  sensitivity  and  specificity  increase 
substantially.  Sensitivity  increases  from  a  range  of  49-88%  for 
either  modality  alone  to  a  range  of  71-89%  with  the  two 
modalities  combined,  while  specificity  increases  from  a  range 
of  57-84%  to  a  range  of  61-91%  (16-18). 

Angiogenesis,  the  formation  of  new  blood  vessels,  is  often 
associated  with  malignancy  (19,20).  A  number  of  studies  have 
reported  an  increased  microvessel  density  (MVD)  in  PCa, 
potentially  enabling  the  detection  or  localization  of  tumors 
through  imaging  techniques  sensitive  to  these  characteristics 
(21-24).  DCE-MRI  of  PCa  is  often  found  to  enhance  more  quickly, 
to  a  greater  degree,  and  to  show  more  washout  than  the  benign 
PZ  (22,25-28).  Malignant  tissue  differs  from  benign  tissue  with  re¬ 
spect  to  microvessel  density,  blood  flow,  vascular  morphology  and 
permeability,  and  flow  dynamics  (29).  A  growing  body  of  literature 
suggests  that  DCE-MRI  may  significantly  improve  cancer  detection, 
tissue  characterization,  localization,  and  staging  (30-33).  According 
to  a  study  by  Kim  et  at.,  DCE-MRI  has  been  found  to  improve  accu¬ 
racy,  sensitivity,  and  specificity  by  26%,  31%,  and  22%  respectively 
relative  to  1.5T  T2W  imaging  alone  (34).  Another  study  by  the 
same  group  at  3T  also  confirmed  the  superiority  of  DCE-MRI  over 
T2W  imaging  (35).  However,  reports  of  correlation  between  DCE- 
MRI  and  specific  tissue  properties  such  as  GS  and  MVD  have  been 
mixed,  with  both  significant  (24,36,37)  and  non-significant  (26,38) 
findings.  Moreover,  mechanisms  governing  both  qualitative  and 
quantitative  changes  are  not  yet  fully  understood. 


MR  SPECTROSCOPIC  IMAGING  OF  PROSTATE 
USING  1 D  SPECTROSCOPY 

MRS  is  a  non-invasive  and  powerful  biochemical  technique,  which 
can  be  performed  after  the  MRI  protocol.  Four  metabolites,  namely 
citrate  (Cit),  creatine  (Cr),  spermine  (Spm),  and  choline  (Ch)  are  the 
ones  commonly  detected  in  PCa  (39,40).  The  resonances  of  the 
above  mentioned  metabolites  occur  at  distinct  frequencies 
(approximately  2.6  ppm,  3.03  ppm,  3.1  ppm,  and  3.2  ppm,  respec¬ 
tively).  In  healthy  prostate  tissues,  Cit  is  secreted  by  the  epithelial 
cells  of  the  prostate  in  large  amounts  along  with  high  levels  of  zinc, 
which  inhibit  the  oxidation  of  Cit  in  the  Krebs  cycle.  In  the  presence 
of  cancer,  the  Cit  level  is  dramatically  diminished  due  to  significant 
reduction  of  zinc  in  the  cancerous  epithelial  tissue.  Concurrently, 
the  Ch  level  is  elevated  due  to  increased  cell  membrane  turnover 
in  the  proliferating  malignant  tissues.  Magnetic  resonance 
spectroscopic  imaging  (MRSI)  is  promising  as  a  valuable  technique 
for  evaluating  the  extent  and  aggressiveness  of  primary  and 
recurrent  PCa  (10,11).  Instead  of  single  voxels,  MRSI  provides 
spectra  from  three-dimensional  (3D)  spatial  arrays  of  contiguous 
volumes  mapping  the  entire  prostate.  MRI  and  MRSI  both  are  used 
for  detailed  anatomic  and  metabolic  evaluations  of  the  prostate. 
The  (Ch  +  Cr)/Cit  ratio  is  usually  used  as  a  criteria  for  the  diagnosis 
of  PCa  in  the  3D  MRSI  studies. 

A  major  limitation  of  the  MRSI  using  ID  spectra  is  its  severe 
spectral  overlap  because  of  limited  spectral  dispersion  at 
clinically  used  static  magnetic  field  strengths  (B0<  3T).  Low 
spectral  dispersion  causes  a  large  number  of  metabolites  to 
overlap  within  a  small  range  of  'H  spectra  (0-5  ppm)  (41). 
Consequently,  it  makes  quantification  of  metabolites  very 
challenging.  Further,  water  and  lipid  suppression  using 
Mescher-Garwood  (MEGA)  RF  pulses  are  combined  with  the 
localizing  point  resolved  spectroscopy  (PRESS)  or  stimulated 
echo  acquisition  mode  (STEAM)  sequence.  This  has  resulted  in 
a  long  TE-based  MRSI  acquisition  since  the  1990s  (10,11,42). 
Spectral  editing  techniques  (43-47)  are  generally  optimized  for 
detecting  one  specific  metabolite,  and  hence  may  not  be 
optimal  for  simultaneously  detecting  a  large  number  of 
metabolites  in  clinical  practice. 

SINGLE  VOXEL  BASED  2D  MR  SPECTROSCOPY 
OF  PROSTATE  CANCER 

Different  versions  of  multidimensional  techniques  have  been 
reported  to  overcome  the  signal  overlap  and  to  detect  several 
metabolites  more  unambiguously  than  conventional  MRSI 
techniques  including  spectral  editing  (46-56).  Multidimensional 
MRS  enables  detection  of  a  larger  number  of  resonances  from 
multiple  metabolites  with  improved  spectral  resolution  than 
one-dimensional  (ID)  MRS  and  a  more  efficient  and  accurate 
identification  and  quantification  of  metabolites.  There  are 
several  books  explaining  various  theoretical  aspects  of  2D  MRS 
(57-59).  Here,  we  present  a  brief  theoretical  explanation  to 
describe  the  acquisition  and  post-processing  of  2D  MRS. 

Theory 

In  the  recording  of  single-voxel-based  1 D  spectroscopy,  once  the 
boundaries  of  the  voxel  are  spatially  localized  with  slice-selective 
RF  pulses,  differences  in  frequencies  of  'H  for  different 
metabolites  are  attributed  to  chemical  shift  or  shielding 
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experienced  by  the  proton  in  the  particular  electronic  environment 
and  indirect  spin-spin  coupling  (J-coupling)  communicated  through 
covalent  bonds.  Chemical  shift  results  in  distribution  of  signal  over 
several  peaks  originating  from  different  magnetically  non-equivalent 
proton  groups.  J-coupling  splits  the  amplitude  of  each  group  further, 
and  distributes  it  over  multiplets  arising  from  indirect  coupling  of 
magnetically  non-equivalent  protons  mediated  by  electron-nuclear 
and  electron-electron  interactions.  Overlap  of  these  multiplets 
makes  quantitation  of  different  metabolites  difficult. 

2D  L-COSY  and  JPRESS 

The  idea  behind  multidimensional  spectroscopy  is  to  pry  open 
these  interactions  and  observe  their  relationships  in  action.  The 
sequences  are  designed  to  see  how  a  spin  group  modifies  its 
response  when  another  is  perturbed  in  some  way.  This  brings  to 
the  foreground  the  correlation  between  interacting  spin  groups, 
which  is  the  basis  of  correlated  spectroscopy  (COSY)  (52,57).  The 
three-dimensional  (3D)  spatially  resolved  analog  of  2D  COSY  was 
named  "L-COSY"  (52).  In  this  scheme,  a  preparation  period  initiates 
a  single  slice-selective  90°  RF  pulse  or  a  sandwich  of  slice-selective 
RF  pulses  in  two  orthogonal  planes  (90°-/)-180°-/l)  during  which 
the  (longitudinal)  equilibrium  magnetization  of  nuclear  spins  is 
transferred  to  transverse  magnetization  or  Flahn  spin-echo.  This 
is  followed  by  encoding  the  second  spectral  dimension  with  a 
variable  time  period,  fq,  meaning  that,  during  a  series  of  repeat 
experiments,  fq  takes  on  a  different  set  of  values.  This  is  usually 
termed  the  evolution  period,  which  is  then  followed  by  a  mixing 
period  containing  a  slice-selective  90°  RF  pulse  in  the  third 
orthogonal  plane  in  2D  L-COSY.  During  this  period,  there  is  a 
coherence  transfer  between  J-coupled  spins.  In  the  localized  2D 
J-resolved  spectroscopy  (JPRESS),  a  slice-refocusing  180°  RF  pulse 
replaces  the  slice-selective  90°  RF  pulse  in  the  L-COSY  sequence 
(49,50).  After  this,  the  data  acquisition  or  detection  period,  f2, 
begins  in  both  L-COSY  and  JPRESS,  during  which  the  digitized 
signal  is  recorded  as  a  function  of  f2,  similarly  to  ID  MRS.  This  is 
repeated  several  times,  creating  a  2D  data  matrix,  with  each  row 
representing  a  different  fq.  This  arrayed  signal  acquisition,  s(f2,  fq), 
is  the  basis  of  2D  spectroscopy  and  can  be  extended  to  further 
dimensions  by  combining  with  two  or  three  spatially  encoding 
gradients,  thereby  enabling  recording  of  multivoxel  2D  spectra  in 
a  single  slice  or  a  3D  volume.  More  spectroscopic  dimensions  will 
come  into  play  when  considering  multinuclear  MRS.  To 
understand  the  nature  of  the  interactions  between  spins  during 
evolution,  mixing,  and  detection  periods,  and  how  these  events 
modulate  the  amplitude,  frequency,  phase,  and  full-width  at  half 
maximum  (FWPIM)  of  the  signal,  we  need  to  take  a  closer  look  at 
the  J-coupled  spin-pair  system. 

A  weakly  (or  strongly)  J-coupled  spin  pair  system  is  defined  as 
one  in  which  the  chemical  shift  differences  between  magnetically 
non-equivalent  protons  (<5|  —  Ss)  within  the  molecule  are  larger  than 
(or  equal  to/less  than)  the  J-coupling  between  the  spins  separated 
by  covalent  bonds.  It  can  be  shown  that  a  general  Flamiltonian  for 
J-coupled  spins  is  given  by  (58) 

H  =  IjCOjtjz  +  Ij<k2Kj)kljlk  [1  ] 

where  a  set  of  spins  (/)  with  different  chemical  shifts  experiences 
mutual  interactions,  and  to  and  J  are  expressed  in  rad/s.  For  a  weakly 
coupled  two-spin  system  [IS,  two  spin-1/?  particles)  with  four 
Zeeman  product  states  in  superposition,  there  are  four  stationary 
states  or  energy  levels  that  represent  the  fractional  population  of 
coupled  spins  at  equilibrium.  The  allowed  transitions  in  a  J-coupled 


spin  pair  system  are  similar  to  that  in  an  isolated  spin  system. 
Flowever,  the  total  magnetic  quantum  number,  m,  of  only  one  spin 
can  change  by  ±1 .  Flence,  there  will  be  four  resonances  containing 
two  doublets  for  the  IS  protons  with  one  for  the  /  and  the  other  for 
the  S  spin. 

2D  L-COSY 

The  2D  signal  from  spin  /  acquired  along  the  detection  dimen¬ 
sion  (f2)  is  given  by  (52) 

s(tq,f2)  =  Tr[(/X)cr]exp^— k»2®t2^exp(— fq/r2)exp(— t2/T2)[1  -  exp(-TR/Tq)] 

[2] 

where 

a«  —  0.5  cos(2ffJ4)[/ycos(cuq^fq)cos(ffJfq)  +  /xsin(roq  Wfq)cos(zrJfq) 
+2/z5)<cos(coq(,lfq)sin(^Jfq)  —  2/zSysin(a)qWfq)sin(;rJfq)] 

[3] 

+0.5  sin(2ffJ/l)[/ycos(coqWfq)sin(^Jfq)  + /xsin(coqWfq)sin(^Jfq) 

— 2/zSxcos(raq('lfq)cos(^Jfq)  +  2/zSysin(cuq  Wfq)cos(^Jfq)] 

Tr,  Tq,  and  T2  represent  repetition  time  and  longitudinal  and 
transverse  relaxation  times,  respectively.  It  is  also  evident  from 
Equation  [3]  that  the  coherence  transfer  (CT)  from  spin  /  to  S  is 
characterized  by  two-spin  operators,  2 lzSx  and  2 lzSy.  A  similar 
equation  can  be  derived  for  5  spin  resulting  in  a  coherence  trans¬ 
fer  to  /  spin. 

Apodization  filters  for  2D  L-COSY 

As  shown  in  Equation  [3],  the  first  maximum  of  the  CT  echo  will 
occur  at  fq  =  1/2  J  for  the  two-spin  system  [IS)  and  the  position  of 
the  maximum  for  lactate  (l3S,  three  methyl  and  one  methine 
protons)  will  be  at  fq  =n/AJ  as  described  previously  by  Ziegler 
et  al.  (60).  In  a  strongly  coupled  AB  spin  system  such  as  citrate 
in  the  prostate,  the  coherence  transfer  function  becomes  more 
complicated,  since  the  chemical  shift  difference  is  equal  to  or 
less  than  J  (50,54,61).  It  is  also  evident  in  Equations  [2]  and  [3] 
that  the  2D  diagonal  peak  intensities  follow  a  cosine 
dependence  and  time  domain  cross-peak  amplitudes  increase 
from  zero  at  the  beginning  to  a  maximum  at  1/2 /with  the  signal 
decay  according  to  T2*.  As  discussed  by  Ernst  et  al.,  for  small  fq 
and  f2  values,  the  contribution  to  the  2D  cross  peak  volume  is 
insignificant  (57).  Flence,  it  is  advisable  to  weight  the  time- 
domain  signal,  as  shown  in  Equation  [2],  by  a  weighting  function 
that  will  deemphasize  the  signal  for  small  f2  and  fq  values  (57). 
While  post-processing  a  2D  L-COSY  spectrum,  optimal  matching 
filters  such  as  a  sine-bell  or  skewed  squared  sine-bell  can  be 
used  along  both  dimensions  for  better  sensitivity  of  2D  cross 
peaks.  A  sine-bell  (SB)  filter  can  be  defined  as 

SBn  =  sin{(*/n)[l  +  (n-  1)f„/(f/max  -  t,min)] }  [4] 

where  n  being  the  shift  parameter  is  a  positive  integer  and  /  runs 
over  the  two  dimensions,  1  and  2.  For  large  n,  Equation  [4] 
reduces  to 

SB  =  sin{7rfw/(f,-max  -  f,min)}  [5] 

As  shown  in  Equation  [5],  the  symmetric  unshifted  SB 
function  has  a  maximum  value  of  1  at  fw  =  (f,max  —  f,mm)/2. 
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As  the  maximum  value  can  be  reached  at  an  earlier  time,  the 
shifted  SB  or  its  square  (SBn2)  may  offer  more  flexibility  than 
the  unshifted  ones. 

As  described  by  Delikatny  et  al.  (62),  there  are  three  major 
advantages  with  using  the  SB  filter.  (1)  As  it  begins  with  zero  value, 
unlike  an  exponential  filter,  it  can  emphasize  cross  peaks  relative  to 
2D  diagonal  peaks  that  are  cosine  dependent.  (2)  It  removes  the 
broad  wings  (dispersive  components)  from  2D  magnitude 
lineshapes.  (3)  At  the  end  of  the  time  domain,  the  trailing  edge  of 
the  SB  function  the  window  function  goes  smoothly  to  zero  and 
truncation  errors  due  to  apodization  are  minimized. 

Shown  in  Figure  1(A)  is  a  2D  L-COSY  spectrum  recorded  in  the 
peripheral  zone  of  a  28  y.o.  healthy  prostate  using  a  1.5T  whole 
body  MRI/MRS  scanner  (GE  Medical  Systems,  Waukesha,  Wl)  with 
an  endorectal  coil  (MEDRAD,  Pittsburgh,  PA)  combined  with  a 
pelvic  phased-array  coil  for  signal  reception.  A  body  RF  coil  was 
used  for  transmitting  the  RF  pulses.  A  2x2x1  cm3  voxel  was 
placed  on  the  lower  left  peripheral  zone  using  an  axial  fast 
spin-echo  MRI.  The  2D  L-COSY  spectrum  was  recorded  using 
the  following  parameters:  TR  =  2s,  minimal  7"E  =  30ms,  45  fq 
increments  and  16  averages.  The  total  acquisition  time  was 
24  min.  The  raw  data  was  acquired  using  1024  complex  points 
and  a  spectral  window  of  2500  Hz  along  the  detected  direct  (f2) 
dimension.  The  incremental  period  (Afq)  was  1.6  ms  to  yield  a 
spectral  window  of  625  Hz  along  the  second  indirect  dimension 
(fq).  The  spectral  2D  raw  files  were  processed  using  the  FELIX 
software  package  (Felix  NMR  Inc.,  San  Diego,  CA).  The  2D  MRS 
array  was  apodized  using  the  skewed  squared  sine-bell  filters 
(62)  along  the  two  axes  and  before  zero-filling  to  2048x256. 
After  double  Fourier  transformation,  the  2D  L-COSY  spectra  were 
reconstructed  in  the  magnitude  mode  and  displayed  as  contour 
plots.  For  comparison,  shown  in  Figure  1 B  is  the  same  2D  L-COSY 
data  processed  using  the  conventional  exponential  filter  func¬ 
tions  (3  Hz),  where  the  emphasis  of  discriminating  2D  cross  peaks 
from  the  diagonal  peaks  is  completely  sacrificed  compared  with 
the  skewed  squared  sine-bell  filters. 


2D  JPRESS 

In  the  half-echo  sampled  2D  JPRESS,  two  increments  (fq/2)  are 
used  before  and  after  the  last  1 80°  RF  pulse  (49,50).  The  2D  signal 
for  spin  /  acquired  along  the  detected  dimension  (f2)  is  the  same 
as  shown  in  Equation  [2]  with  a  redefined  a  and  coq.  As  there  is 
no  mixing  period  here,  the  2D  JPRESS  spectrum  contains  the 
same  number  of  peaks  as  the  conventional  ID  MRS  using  the 
PRESS  sequence.  The  second  spin-echo  during  which  the  second 
spectral  dimension  is  encoded  will  refocus  the  chemical  shift  but 
not  in  bilinear  interaction,  J,  similar  to  the  first  Hahn  echo.  The 
frequency  of  each  peak  along  the  fq  dimension  will  be 
dependent  on  J  only,  where  the  frequency  along  the  f2 
dimension  will  be  (chemical  shift±1/2J).  A  more  detailed 
theoretical  explanation  is  discussed  by  Thrippleton  et  al.  (63). 

2D  JPRESS  of  a  healthy  human  prostate 

Figure  2(B)  shows  a  2  ml  voxel  2D  JPRESS  spectrum  of  the  same 
28  y.o.  healthy  prostate  as  shown  in  Figure  1.  The  voxel  location 
of  the  2D  JPRESS  spectrum  is  shown  in  Figure  2(A).  The  2D  peaks 
due  to  Cit  were  located  along  Fq  =  ±1 .6  Hz,  ±7.9  Hz,  and  ±1 7.5  Hz. 
In  addition,  the  triplet  nicely  resolved  about  F2  =  3.1  ppm  along 
Fq=0Hz  and  ±7.8  Hz  was  identified  as  Spm.  The  presence  of 
strong  Cit  peaks  was  consistent  with  the  well-known  fact  of  its 
high  abundance  in  healthy  prostate  (39,40).  Moreover,  there 
were  2D  peaks  along  Fq  =  0  Hz  due  to  Cr  and  Ch  at  F2  =  3.0  ppm 
and  3.2  ppm,  respectively.  The  J-resolved  peaks  due  to  Ch 
methylene  protons  were  not  observable  due  to  the  smaller  voxel 
size  and  reduced  number  of  averages. 

Apodization  filters  for  2D  JPRESS 

The  J-coupled  multiplets  are  better  resolved  along  the  fq 
dimension  than  the  detected  f2  dimension  as  well  the  ID  MRS 
counterpart,  since  any  defocusing  linear  B0  interactions  including 
the  static  field  inhomogeneities  during  the  first  half  of  fq  are 


Figure  1.  A  1.5  T  2D  L-COSY  spectrum  recorded  from  a  healthy  prostate  of  a  28  y.o.  male  using  two  different  apodization  filters:  (A)  a  skewed  squared 
sine-bell  and  (B)  an  exponential  filter  using  a  3  Hz  line  broadening.  The  2D  L-COSY  spectrum  was  recorded  using  the  following  parameters.  A 
2x2x1  cm3  voxel  was  placed  on  the  lower  left  peripheral  zone,  TR  =  2s,  minimal  TE  =  30  ms,  45  fq  increments,  incremental  period  (Afq)  of  1.6  ms  to  yield 
a  spectral  window  of  625  Hz  along  the  F,  dimension,  16  averages  per  Afq,  1024  complex  points  and  a  spectral  window  of  2500  Hz  for  the  F2  dimension. 
The  total  acquisition  time  was  24  min.  An  endorectal  coil  combined  with  a  pelvic  phased-array  coil  for  signal  "receive"  and  a  body  rf  coil  for  "transmit" 
were  used.  The  raw  spectral  files  were  processed  using  the  FELIX  software  package  (Felix  NMR  Inc.,  San  Diego,  CA).  After  zero-filling  to  2048x256, 
double  Fourier  transformation  was  performed  and  the  2D  L-COSY  spectra  were  reconstructed  in  the  magnitude  mode  and  displayed  as  contour  plots. 
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Figure  2.  (A)  An  axial  fast  spin-echo  MRI  slice  showing  the  2D  MRS  voxel;  the  MRI  acquisition  parameters  were  as  follows:  4  mm  slice,  TR  =  2.5  s,  TE  = 
84ms,  FOV=  14-24cm,  acquisition  matrix  256x  192,  and  4  number  of  excitations  (NEX),  resulting  in  an  acquisition  time  of  4min.  (B)  A  half-echo 
sampled  1.5T  2D  JPRESS  spectrum  recorded  in  the  peripheral  zone  of  the  same  28  y.o.  healthy  prostate  as  used  for  Figure  1  using  the  following 
parameters:  2x2x1  cm3  voxel,  TR=  2  s,  minimal  TE  =  30ms,  45  q  increments,  incremental  period  (Afn)  of  10  ms  (5  ms  before  and  5  ms  after  the  last 
180°  RE  pulse)  to  yield  a  spectral  window  of  ±50Flz  along  the  F,  dimension,  16  averages  per  Af,,  1024  complex  points  and  a  spectral  window  of 
2500  Flz  for  the  F2  dimension.  The  total  acquisition  time  was  24  min. 


refocused  during  the  second  half,  resulting  in  a  net  zero 
dependence  on  the  B0  static  field  inhomogeneities  and  other 
linear  interactions.  Even  though  this  is  a  major  advantage,  the 
phase-modulated  time  domain  datasets  are  transformed  into 
phase-twisted  2D  peaks  after  the  double  fast  Fourier 
transformation  (FFT)  of  the  2D  JPRESS  raw  data.  Flence,  the 
above  mentioned  squared  or  simple  sine-bell  filter  functions 
can  be  used  prior  to  the  double  FFT. 

Strong  coupling  effects  in  2D  JPRESS 

It  was  demonstrated  earlier  that  2D  JPRESS  spectra  of  brain  and 
prostate  metabolites  show  more  cross  peaks  than  those  of 
weakly  coupled  ones  (50,54).  Shown  in  Figure  3  is  a  simulated 
hard  pulse  version  of  the  2D  JPRESS  (90°-180°-180°)  spectrum 
of  Cit  using  the  GAMMA  simulation  library  (64).  Cit  has  two 
equivalent  methylene  groups.  Each  proton  pair  forms  a  strongly 
J-coupled  AB  spin  system  (50,54),  resulting  in  eight  J-resolved  2D 
peaks  anti-symmetric  about  F2  =  2.65  ppm,  as  evident  in  the 
experimental  and  simulated  2D  spectra.  The  2D  peaks  located 
along  Fq  =  ±1 .6  Flz,  ±7.8  Flz,  and  ±1 6.4  Flz  were  in  agreement  with 
a  previous  report  (54).  Besides,  the  projected  1 D  spectra  onto  Fq 
and  F2  axes  are  also  shown. 

2D  JPRESS  of  PCa  patients  with  two  different  Gleason  scores  (GS) 

Figure  4  shows  2D  JPRESS  spectra  recorded  in  a  59  y.o.  BPFH  and 
a  50  y.o.  PCa  patient  using  a  Siemens  1.5T  Avanto  MRI  scanner. 
The  metabolites  such  as  Cit,  Ch,  Cr,  and  Spm  were  identified 
and  the  2D  multiplet  patterns  of  Cit  and  Spm  in  JPRESS  spectra 
were  detected  with  a  reasonable  resolution.  Flowever,  the 
limited  spectral  resolution  along  the  second  axis  (Fq)  resulted  in 
an  overcrowded  2D  JPRESS  spectrum.  This  complex  spectrum 
is  due  to  the  strong  coupling  effects  of  Cit  as  presented  in 
Figure  3.  Nagarajan  et  at.  (56)  have  shown  a  decrease  of  Spm 
in  PCa  patients  of  higher  GS  (4  +  3)  compared  with  lower  GS 


IL 


F2  (ppm) 

Figure  3.  A  simulated  1.5T  2D  JPRESS  spectrum  of  citrate  using 
7=1 5.6 Hz  and  <5  =  9.6 Flz  at  1.5T  (TR  =  2s,  TE  =  30ms,  256 F,  points,  and 
4096  F2  points). 


(3  +  4)  using  2D  J-resolved  spectroscopy.  FHigh  levels  of  Spm  are 
found  in  the  healthy  prostatic  ducts  (65-67)  and  the  observed 
variations  of  Spm  in  cancer  may  be  due  to  the  loss  of  ductal 
morphology  or  a  reduction  in  the  secretion  of  polyamines  (68). 
The  Spm  itself,  which  plays  a  role  in  cell  proliferation  and 
differentiation,  may  provide  additional  information  for  early 
diagnosis  and  prognosis  predicting  tumor  progression.  Flowever, 
the  separation  of  Ch  from  Spm  is  inherently  difficult  because  of 
the  proximity  of  the  peaks  and  the  limited  spectral  resolution 
available  within  a  reasonable  imaging  time.  A  distinct  integration 
region  cannot  be  assigned  for  Spm.  Thus,  polyamine  level  could 
not  be  integrated  and  quantified  unambiguously  using  operator 
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Figure  4.  T2-weighted  axial  MRI  slice  images  showing  locations  of  2D  JPRESS  spectra  recorded  in  the  peripheral  zones  of  (A,  C)  a  59  y.o.  BPH  patient 
and  (B,  D)  a  50  y.o.  prostate  cancer  patient.  The  half-echo  sampled  2D  JPRESS  spectra  were  recorded  using  the  following  parameters:  a  1.5T  MRI 
scanner,  TR  =  2s,  minimal  TE  =  30ms,  64f,  increments,  incremental  period  (At,)  of  10ms  (5  ms  before  and  5ms  after  the  last  180°  RF  pulse)  to  yield  a 
spectral  window  of  ±50  Hz  along  the  F,  dimension,  eight  averages  per  Afq,  2048  complex  points  and  a  spectral  window  of  2000  Hz  for  the  F2  dimension. 
The  total  acquisition  time  was  16  min.  The  raw  spectral  files  were  processed  using  the  FELIX  software  package  (Felix  NMR  Inc.,  San  Diego,  CA).  After 
zero-filling  to  2048x256,  double  Fourier  transformation  was  performed  and  the  2D  JPRESS  spectra  were  reconstructed  in  the  magnitude  mode  and 
displayed  as  contour  plots. 


defined  peak  volumes  in  the  frequency  domain.  It  can  be  only 
assessed  qualitatively  and  reported  relative  to  the  choline  peak.  In 
the  prostate,  Spm  has  three  multiplets  at  1.8  ppm,  2.1  ppm,  and 
3.05-3.15  ppm.  The  3.05-3.1 5  ppm  multiplet  occurs  between  the 
Cr  and  Ch  singlets  and  cannot  be  entirely  resolved  from  them. 
Although  this  ratio  has  been  traditionally  called  the  (Ch+Cr)/Cit 
ratio,  it  actually  includes  Spm  also.  This  presents  a  problem  when 
Ch  is  elevated  in  cancer  and  Spm  is  diminished.  This  hinders 
accurate  quantification  of  both  metabolites.  Application  of  localized 
2D  JPRESS  was  also  reported  in  prostate  in  vivo  studies  by  another 
group  (69). 

Maximum-echo  sampled  2D  JPRESS 

In  the  half-echo  sampled  2D  JPRESS,  if  the  incremental  period 
(A fq)  typically  used  is  10  ms,  the  resulting  total  echo  time 
duration  will  be  630  ms  to  achieve  a  spectral  window  of  100  Hz 
along  F,..  This  would  lead  to  severe  T2  attenuation  of  2D  peaks, 
resulting  in  a  significant  signal  loss.  It  is  known  that  T2  value  of 
several  metabolites  are  shorter  at  3T  than  1.5T  (70,71),  which 


would  lead  to  further  attenuation  of  J-resolved  2D  peaks  at 
higher  field  strength.  As  reported  by  Schulte  et  al.  recently,  a 
maximum-echo  based  JPRESS  can  sample  the  echo  signals 
starting  immediately  after  the  final  crusher  gradient  of  the  last 
180°  pulse  (72).  This  acquisition  scheme  has  several  advantages 
over  the  half-echo  sampling  data  presented  in  the  previous 
section,  where  the  acquisition  starts  at  the  echo  top. 

A  common  definition  of  the  sensitivity  is  the  signal-to-noise 
(SNR)  ratio  per  unit  time  for  identical  experimental  durations 
and  a  constant  acquisition  time  window;  the  noise  will  remain 
the  same  and  it  is  sufficient  to  compare  the  signal.  The  signal 
is  commonly  defined  as  the  peak  height  in  the  frequency 
domain,  which  is  equivalent  to  the  integral  of  the  time  domain 
signal.  For  single  resonances,  it  suffices  to  integrate  over  the 
exponential  damping  curve  of  the  echo  along  the  times  fq  and  t2. 

As  described  by  Schulte  et  al.,  the  efficiency  of  2D  experiments 
can  be  compared  with  1 D  PRESS  by  integrating  over  the  indirect 
dimension  (fq).  The  maximum-echo  sampled  JPRESS  signal  can 
be  divided  into  two  halves  (72).  The  damping  curve  in  the 
traditional  half-echo  sampling  JPRESS  is  equivalent  to  the  right 
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half  as  shown  in  the  following  equation: 

|Tslrexp(-f2/r2*)  exp(— fi/T2)  dt2dfi 
=  r2*T[1  -exp(-Tsl/r2)]  [6] 

where  Ts\  represents  the  total  sampling  time  along  fq. 

Adding  both  halves  together,  the  total  signal  of  the  maximum- 
echo  sampled  JPRESS  can  be  written  as 

2T2*t2[i  -  exp(-rs,/r2)]  -  r2*tx[i  -  exP(-rs,/rx)]  [7] 

where 

fx  =  2  T2T2*/(T2  +  2  T2*)  [8] 

In  contrast,  the  equation  for  the  ID  PRESS  sequence  with  the 
shortest  TE  equivalent  to  the  2D  JPRESS  at  fq  =0  can  be  written 
as  the  following: 

fsrexp(-t2/T2*)  df2dfq  =  T2* TsI  [9] 

The  ratio  of  maximum-echo  sampled  JPRESS  to  PRESS  is  given 
by 

2t2/ts, [i  -  exP(-rs,/r2)]  -  fx/rs,  [i  -  exp(-rsl/fx)]  [io] 
and  half-echo  JPRESS  to  PRESS  (Equation  [9])  by 

r2/rs[i  -  exp(-rsl/r2)]  [ii] 

As  presented  by  Schulte  ef  at.,  the  maximum-echo  sampled 
JPRESS  scheme  has  optimal  sensitivity  in  2D  experiments  and  it 
has  the  same  sensitivity  as  the  ID  PRESS  at  long  72>200ms, 
which  is  the  case  for  singlets  from  Cr  and  trimethyl  protons  of 
Ch  (72). 

Shown  in  Figure  5  is  maximum-echo  sampled  2D  JPRESS  data 
acquired  in  a  27  y.o.  healthy  male  before  (A)  and  after  (B)  a  phase 
rotation  using  the  following  parameters:  a  voxel  of  2  x  2  x  2  cm3, 
fR/TEmin  =  2s/30ms,  2048  complex  points  for  the  detected  f2 
dimension  sampling  a  spectral  width  of  2000  Hz,  100  Afq 
increments  before  the  last  1 80°  slice-selective  RF  pulse  with  each 
Afq  =  1  ms,  eight  averages  per  Af1(  a  3T  MRI  scanner,  16  channel 


body  matrix  coil  for  "receive"  and  a  quadrature  body  "transmit" 
coil.  This  is  in  agreement  with  what  was  shown  by  Lange  et  at. 
(69). 

ECHO-PLANAR  CORRELATED  AND 
J-RESOLVED  SPECTROSCOPIC  IMAGING 

The  single-voxel  (SV)-  based  ID  and  2D  MRS  studies  suffer  from 
limited  spatial  coverage  due  to  recording  of  one  voxel  per 
measurement.  In  contrast,  the  2D  or  3D  MRSI  technique  would 
facilitate  adequate  spatial  coverage  in  a  single  recording. 
Flowever,  a  major  limitation  of  the  conventional  2D/3D  MRSI 
technique  stems  from  using  incremented  phase  encoding  for 
two  or  three  directions  to  traverse  the  /(-space  and  the  total  scan 
time  required  for  the  acquisition  of  a  high-resolution  3D  MRSI 
data  may  be  prohibitively  long  for  clinical  exams.  Another 
drawback  is  that  only  four  major  metabolites  (Cit,  Cr,  Spm,  and 
Ch)  have  been  detected  by  the  3D  MRSI  sequence  due  to  long 
TE  for  optimal  suppression  of  water  and  lipids.  The  acquisition 
of  fully  phase-encoded  MRSI  can  be  greatly  shortened  by  using 
echo-planar  spectroscopic  imaging  (EPSI),  a  method  originally 
proposed  by  Mansfield  (73),  in  which  a  time  varying  readout 
gradient  encodes  spatial  and  spectral  dimensions  during  a  single 
readout.  Much  effort  has  been  devoted  to  implementing  EPSI 
(74-78),  the  results  of  which  have  shown  sufficient  SNR.  In  the 
last  few  years,  two  different  fast  MRSI  sequences  have  been 
demonstrated  with  potential  applications  in  PCa.  First,  the 
flyback  echo-planar  read-out  trajectories  were  incorporated  by 
Chen  et  at.  into  the  PRESS  sequence  using  composite  RF  pulses 
with  high-quality  MRSI  data  recorded  in  nine  PCa  patients  (79). 
Second,  high-resolution  EPSI  was  implemented  on  a  4.7 T  MRI 
scanner  with  the  sequence  tested  in  animal  models  of  PCa  (80). 

To  overcome  a  major  limitation  of  the  SV  localized  2D  MRS, 
phase-encoding  gradients  can  be  combined  with  L-COSY  and 
JPRESS  to  record  multivoxel  2D  spectra  in  human  tissues. 
Flowever,  the  total  duration  will  be  impractically  long  due  to  four 
different  increments  for  encoding:  one  for  the  second  spectral 
dimension  and  three  for  the  three  spatial  dimensions.  With  the 
above  mentioned  progress  in  accelerating  one  spatial  and 
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Figure  5.  A  maximum-echo  sampled  2D  JPRESS  spectrum  recorded  in  the  peripheral  zone  of  A  27  y.o.  healthy  prostate  using  a  3T  MRI  scanner  after 
(A)  a  double  FFT  of  the  2D  raw  matrix  and  (B)  a  double  FFT  after  incorporating  the  phase  rotation  of  the  raw  matrix  to  impose  the  evolution  after  the  last 
180°  slice-selective  RF  pulse. 


NMR  Biomed.  (2013) 


Copyright  ©  2013  John  Wiley  &  Sons,  Ltd. 


wileyonlinelibrary.com/journal/nbm 


NMR 

IN  BIOMEDICINE 


M.  A.  THOMAS  ET  AL. 


spectral  coverage  using  the  3D  EPSI,  the  2D  L-COSY  and  JPRESS 
sequences  were  recently  modified  with  an  EPSI  type  of  readout  to 
yield  2D  spectra  from  multiple  voxels  in  a  single  experiment,  called 
echo-planar  correlated  spectroscopic  imaging  (EP-COSI)  (81)  and 
echo-planar  J-resolved  spectroscopic  imaging  (EP-JRESI)  (82), 
respectively.  The  EP-COSI  and  EP-JRESI  sequences  use  a  bipolar 
spatial  read-out  train  facilitating  simultaneous  spatial  and  spectral 
encoding,  and  the  conventional  phase  and  spectral  encodings  for 
the  remaining  spatial  and  indirect  spectral  dimensions.  Multiple  2D 
COSY  and  2D  J-resolved  spectra  were  recorded  over  the  spatially 
resolved  volume  of  interest  (VOI)  localized  by  a  train  of  three  slice-se¬ 
lective  RF  pulses  (90°-1 80°-90°)  and  (90°-180°-180°).  These  two 
techniques  can  be  used  to  detect  and  quantify  metabolites  with  less 
T2  weighting  than  the  earlier  3D  MRSI  sequences  using  conventional 
phase  encoding  because  of  the  shorter  echo  time.  Figure  6  shows 
the  four-dimensional  (2D  spectral  +  2D  spatial)  EP-COSI  data 
recorded  in  a  27  y.o.  healthy  volunteer.  The  voxel  placement  is 
shown  in  Figure  6(A)  and  the  multivoxel  display  of  Cr,  Ch,  and 
Spm  over  the  localized  VOI  in  Figure  6(B).  An  extracted  3  ml  2D  COSY 
spectrum  is  displayed  in  Figure  6(C). 

A  500  ml  prostate  phantom  was  prepared  containing  the 
following  metabolites  at  physiological  concentrations  as  reported 
in  healthy  human  prostate  (69):  Cit  (50  mM),  Cr  (5  mM),  Ch  (1  mM), 
Spm  (6mM),  ml  (10mM),  PCh  (2mM),  taurine  (Tau,  3  mM),  Glu 
(4  mM),  Gin  (2.5  mM)  and  scyllo-inositol  (Scy,  0.8  mM).  Here,  we  show 
that  all  of  the  above  mentioned  metabolites  can  be  detected  using 
fully  encoded  EP-JRESI  data.  Shown  in  Figure  7(A)  is  an  axial  MRI  slice 
image  showing  the  multivoxel  grids  of  MRSI  with  the  yellow 
boundary  of  the  field  of  view  (FOV),  and  the  white  box  representing 
the  volume  of  interest  (VOI)  localized  by  the  PRESS  sequence  which 
is  an  integral  part  of  the  EP-JRESI  sequence  (82).  The  following 
parameters  were  used  for  acquiring  the  4D  EP-JRESI  data:  TR/TE  = 
1500ms/30ms,  16  phase  encodes  (ky),  32  read-out  points  ( kx )  with 
oversampling,  51 2  ±  read-out  trains  resulting  in  512  pairs  of  complex 
spectral  points  in  the  second  spectral  dimension  (f2),  100fq 
increments  for  the  indirect  spectral  dimension  and  one  average 
per  encoding.  A  total  duration  of  40  min  was  necessary  to  acquire 
this  water-suppressed  4D  EP-JRESI  data.  A  non-water-suppressed 
EP-JRESI  data  using  four  averages  with  only  one  fq  increment  was 
used  for  eddy  current  and  phase  correction  of  the  suppressed  data 
(81).  After  apodization  and  Fourier  transformation  of  this  4D  data, 
the  reconstructed  2D  J-resolved  spectra  were  overlaid  on  top  of 


the  16x16  spatial  grids.  An  extracted  2D  J-resolved  spectrum 
(3  ml)  around  the  center  of  the  VOI  is  shown  in  Figure  7(B).  The 
diagonal  peaks  cutting  through  Fq  =0  show  all  singlets  contained 
in  a  1 D  PRESS  spectrum  and  the  2D  cross  (off-diagonal)  peaks  of 
J-coupled  metabolites  such  as  Cit,  Spm,  Glu,  Gin,  and  ml  were  also 
clearly  visible. 

MULTI  ECHO  (ME)  BASED  ECHO-PLANAR 
J-RESOLVED  SPECTROSCOPIC  IMAGING 
(MEEP-JRESI) 

Multi-echo  (ME)  encoding  schemes,  namely  turbo  spin  echo 
(TSE)  and  fast  spin  echo  (FSE),  have  been  shown  to  decrease 
the  overall  scan  time  in  MRI  (83-85).  Similar  approaches 
have  demonstrated  the  applicability  of  ME  techniques  to 
MRSI  (86,87).  One  limiting  factor  for  ME-based  MRSI  is  the 
T2  decay  as  the  signal  is  greatly  diminished  with  each  echo, 
especially  in  living  tissues,  where  T2  relaxation  times  are 
shorter  (70,71,88).  Recently,  a  4D  EP-JRESI  sequence  was 
implemented  combining  two  spectral  dimensions  with  two 
spatial  dimensions  and  incorporating  multi-echo  (ME)  for 
encoding  one  of  the  spatial  dimensions  to  reduce  scan 
times  to  suit  clinical  requirements.  Nagarajan  ef  al.  (89) 
employed  the  MEEP-JRESI  sequence  on  a  3T  MRI/MRS 
scanner  and  evaluated  it  in  three  healthy  prostate 
volunteers  using  the  external  body  matrix  "receive'  coil. 
The  MEEP-JRESI  technique  facilitates  recording  multivoxel 
2D  J-resolved  spectra  in  a  single  recording  using  a  total 
acquisition  time  of  approximately  13  min.  Figure  8(A)  shows 
the  r2-weighted  axial  MRI  displaying  the  VOI  location  for  the 
MEEP-JRESI  recorded  in  a  28  y.o.  healthy  volunteer.  A 
multivoxel  display  of  Cr/Ch/Spm  inside  the  VOI  is  shown  in 
Figure  8(B)  with  each  voxel  resolution  of  2  ml. 

4D  EP-JRESI  AND  EP-COSI:  ACCELERATED 
ACQUISITION  AND  COMPRESSED  SENSING 
RECONSTRUCTION 

Image  acquisition  approaches  have  conventionally  followed  the 
Nyquist-Shannon  sampling  theorem,  where  the  sampling  rate 


(A) 


Figure  6.  (A)  T2-weighted  axial  MRI  of  27  y.o.  healthy  volunteer  with  MRSI  voxel  location  using  a  3T  MRI  scanner  and  (B)  multivoxel  distribution  of  Cr, 
Ch,  and  Spm  peaks  extracted  from  the  EP-COSI  data.  (C)  2D  L-COSY  spectrum  extracted  from  the  EP-COSI  data  of  a  27  y.o.  healthy  male. 
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Figure  7.  (A)  An  axial  slice  MRI  of  a  prostate  phantom  containing  10  metabolites;  (B)  an  extracted  2D  J-resolved  spectrum  (voxel  size  of  3  ml).  A  3T  MRI 
scanner  was  used. 


(A)  (B) 

Figure  8.  (A)  The  T2-weighted  axial  MRI  and  (B)  multivoxel  distribution  of  Cr,  Ch,  and  Spm  from  a  28  y.o.  healthy  prostate  from  MEEP-JRESI.  A  3  T  MRI 
scanner  was  used. 


should  be  at  least  twice  the  maximum  frequency  present  in  the 
signal  of  interest.  However,  the  theory  of  compressive  sensing 
(CS)  describes  a  combined  sampling  and  reconstruction 
paradigm  which  states  that  certain  images  and  signals  can  be 
recovered  from  an  acquisition  that  uses  fewer  samples  than 
required  by  Nyquist-Shannon  (90,91).  Since  the  annual  meeting 
of  the  International  Society  for  Magnetic  Resonance  in  Medicine 
in  May  2007,  CS  has  met  with  significant  enthusiasm  among  MR 
researchers  (92-96).  Three  criteria  are  critical  for  the  successful 
application  of  CS  (92):  (1)  the  data  should  have  a  sparse 
representation  in  a  particular  transform  domain,  (2)  the  aliasing 
from  undersampling  should  be  incoherent  in  that  transform 
domain  and  (3)  a  nonlinear  reconstruction  should  be  used  to 
enforce  both  the  sparsity  of  the  reconstruction  and  consistency 
with  the  measurements.  MRI  is  well  suited  for  CS,  and  there  are 
significant  benefits  in  imaging  speed  and  reduced  costs,  thereby 
improving  patient  care.  A  major  challenge  in  designing  CS  data 
acquisition  methods  for  MRI  is  in  implementing  NUS  densities 
that  result  in  incoherent  aliasing  while  providing  data  sparsity 
in  a  transform  domain,  such  as  wavelets,  curvelets,  etc  (76). 
Incoherent  aliasing  combined  with  sparsity  in  the  transform 


domain  allows  LI -norm-based  reconstructions  from  NUS  data  to 
be  exact  under  ideal  conditions  and  "approximately"  exact  under 
normal  conditions  (90,91).  Application  of  CS  sampling  and 
reconstruction  has  been  accomplished  in  parallel  imaging  to 
exploit  both  image  sparsity  and  coil  sensitivity  encoding  (97).  As 
described  earlier,  there  are  three  spatial  encodings  and  one 
spectral  encoding  in  a  3D  MRSI  sequence,  and  recently  Vigneron 
and  co-workers  developed  an  undersampling  scheme  along  both 
spatial  encoded  dimensions  to  achieve  suitable  incoherent 
aliasing.  They  have  demonstrated  a  factor  of  two  enhancement 
in  the  spatial  resolution  without  increasing  acquisition  time  or 
decreasing  coverage  (94).  A  further  modified  scheme  was  shown 
to  provide  up  to  a  acceleration  factor  of  7.5  for  hyperpolarized 
MRSI  (98). 

The  total  duration  of  a  fully  encoded  3D  MRSI  can  be 
approximately  2  hours  if  the  following  parameters  are  used:  TR  = 
Is,  12x12x12  for  phase-encoding  along  three  spatial  dimen¬ 
sions  and  four  averages.  Therefore,  average  weighted  and  other 
schemes  have  been  used  to  minimize  the  total  time  (99).  During 
the  last  five  years,  two  different  fast  MRSI  sequences  have  been 
demonstrated  for  PCa  as  mentioned  earlier  (79,80).  There  has 
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been  one  report  of  combining  two  spectral  dimensions  with  2D 
spatial  encoding  applicable  to  PCa  (100).  As  discussed  earlier,  a 
second  spectral  dimension  was  added  to  the  EPSI  technique  in 
the  recently  published  4D  EP-COSI  sequence  (81).  A  total  scan 
time  of  40  min  was  required  for  51 2  x  1 00  complex  points  along 
the  two  spectral  dimensions  with  a  16x16  spatial  resolution 
along  the  two  spatial  dimensions. 

We  demonstrate  here  that  by  imposing  a  NUS  scheme  on  the 
above  presented  fully  sampled  data  shown  in  Figure  7(B)  and 
after  the  CS  reconstruction,  the  quality  of  the  extracted  2D 
MRS  spectra  from  the  undersampled  data  is  comparable  to  that 
of  the  fully  sampled  data.  Therefore,  NUS  was  imposed  on  one 
spectral  and  one  spatial  dimension  and  the  remaining  spectral 
and  spatial  dimensions  were  fully  sampled  by  the  EPI  readout. 
Compared  with  the  fully  encoded  data,  the  25%  NUS  densities 
demonstrate  4x  undersampling,  which  reduced  total  scan 
duration  to  10  min.  Figure  9  shows  the  multivoxel  EP-JRESI  data 
(A)  and  extracted  2D  J-resolved  spectrum  (B)  after  the  CS  recon¬ 
struction,  demonstrating  the  fidelity  of  the  CS  reconstruction. 

The  feasibility  of  recording  a  25%  sparsely  sampled  in  vivo 
EP-JRESI  data  and  successful  CS  reconstruction  of  the  2D 
J-resolved  spectrum  is  demonstrated  here.  A  32  y.o.  healthy 
human  subject  was  investigated  on  the  3T  MRI  scanner 
using  the  quadrature  body  coil  "transmit"  and  external  body 
matrix  "receive"  coil  assembly.  The  following  parameters 
were  used  to  acquire  the  NUS  sampled  EP-JRESI  data:  TR/ 
TE  =  1 .5  s/30  ms,  2  averages,  512  f2,  oversampled  32  kx,  25% 
NUS  along  the  indirect  spectral  (fn )  and  spatial  ky  dimen¬ 
sions.  In  Figure  10(A),  an  axial  MRI  of  the  abdomen  is  shown 
displaying  the  VOI  covering  the  prostate  localized  by  the 
PRESS  sequence  and  the  EP-JRESI  grids.  The  split-Bregman 
iterative  algorithm  (101)  was  used  to  reconstruct  the 
missing  fn  and  ky  data  points  from  the  prospectively 
undersampled  data  matrix.  The  CS-reconstructed  multivoxel 
display  confirmed  the  feasibility  of  detecting  prostate 
metabolites  over  many  voxels.  The  CS-reconstructed  2D 
J-resolved  spectrum  extracted  from  the  central  location  is 
shown  in  Figure  10(B).  The  2D  diagonal  and  cross  peaks  of 
Cit  and  other  metabolites  are  visible  as  reported  in  a  recent 
publication  (102).  The  endorectal  "receive"  coil  is  expected 


to  facilitate  at  least  one  order  of  magnitude  higher  sensitivity  than 
that  of  the  body  matrix  assembly  used  here.  Hence,  significantly 
improved  sensitivity  can  be  achieved  using  the  endorectal  coil. 
Our  recent  findings  on  the  endorectal  coil-based  NUS  EP-JRESI 
acquisition  and  CS  reconstruction  confirming  the  improved 
detectability  of  prostate  metabolites  will  be  published  elsewhere. 

PRIOR-KNOWLEDGE  FITTING  FOR 
QUANTITATION  OF  PROSTATE  METABOLITES 

As  discussed  earlier,  previous  attempts  by  other  researchers  using 
the  single-  and  multivoxel-based  1 D  MRS  approaches  have  reported 
four  prostate  metabolites  only,  such  as  Cit,  Ch,  Spm,  and  Cr  (10,1 1). 
This  was  due  to  the  long  echo  time  (TE>  100ms)  and  limited  ID 
spectral  quantitation  approaches  used.  Recent  investigations  of 
HR-MAS  of  ex  vivo  prostate  cancer  specimens  on  ultra-high  field 
NMR  spectrometers  have  demonstrated  quantitation  of  many  more 
metabolites  such  as  Tau,  ml,  scyllo-inositol  (si),  Glx,  etc.  than  what 
has  previously  been  shown  using  the  ID  MRS  in  vivo  approaches 
so  far  (103-105).  A  recent  report  has  used  prior-knowledge-based 
LC-model  processed  HR-MAS  ex  vivo  data  to  quantify  additional 
metabolites  such  as  Glx  and  glucose  (106).  A  few  years  ago,  Schulte 
et  al.  developed  a  prior-knowledge  fitting  (ProFit)  algorithm  and 
demonstrated  the  feasibility  of  quantitation  of  brain  and  prostate 
metabolites  (69,72,107).  After  the  fitting  process,  the  quality  of  the 
fit  can  individually  be  evaluated  for  each  metabolite  using 
Cramer-Rao  lower  bounds  (CRLB)  (108).  A  statistical  lower  bound 
for  the  achievable  standard  deviation  of  the  estimated  parameters 
is  provided  by  CRLB  which  are  not  dependent  on  the  individual 
concentrations  but  only  on  the  noise  and  the  orthogonality  of  the 
basis  function.  The  architecture  of  the  fitting  process  allows  for 
another  useful  measure  of  the  quality  of  the  fitting  of  the  spectrum 
by  comparing  creatine  3.9  (Cr3.9)  to  creatine  3.0  (Cr3.0)  ratios,  which 
ideally  should  be  1  since  the  number  of  protons  (2,3)  are  already 
taken  into  account  for  Cr3.9  and  Cr3.0,  respectively,  in  their  basis- 
set  creation.  Higher  Cr3.9/Cr3.0  ratios  reflect  poor  spectra  and 
implying  that  the  results  can  be  excluded.  Note  that  in  order  to 
implement  this  control  creatine  3.9  and  creatine  3.0  have  to  be 
implemented  as  different  elements  in  the  basis  sets.  Our  preliminary 
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Figure  9.  (A)  The  CS-reconstructed  2D  J-resolved  multivoxel  spectra  showing  the  Cit  multiplets  processed  from  the  25%  NUS  raw  data  using  a  3  T  MRI 
scanner.  (B)  the  corresponding  2D  J-resolved  spectrum  extracted  from  the  center  voxel. 
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Figure  10.  (A)  An  axial  MRI  slice  of  the  32  y.o.  healthy  male  subject  showing  the  VOI  and  MRSI  grids.  (B)  The  2D  J-resolved  spectrum  extracted  from  a 
2  ml  voxel.  A  3T  MRI  scanner  was  used. 
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Figure  11.  ProFit  fitting  of  2D  JPRESS:  (a)  a  2D  JPRESS  spectrum 
recorded  in  the  peripheral  zone  of  the  27  y.o.  healthy  prostate  using  a 
3  T  MRI  scanner,  (b)  its  fit  and  (c)  the  fit  residue  as  determined  with  ProFit. 


spectra  in  multiple  regions.  However,  the  total  duration  was 
typically  20-40  min  depending  on  the  number  of  incremented 
spectral  (fq)  and  spatial  encoding  (ky)  steps.  Recent  developments 
from  our  group  demonstrate  clearly  that  further  acceleration  can 
be  accomplished  using  the  NUS  schemes  and  total  durations  of 
the  4D  EP-JRESI  and  EP-COSI  sequences  can  be  shortened  to 
approximately  10  min  or  less.  The  nonlinear  CS  reconstruction  is 
necessary  for  processing  the  NUS  4D  data.  Further  work  is 
necessary  to  demonstrate  the  potential  of  multidimensional  MR 
spectroscopic  imaging  using  the  fast  imaging  methods  as  reported 
recently  to  bring  it  into  the  clinic  as  a  robust  diagnostic  metabolite 
imaging  technique.  We  hope  all  these  recent  developments  will 
lead  to  clinical  realization  of  these  novel  MRSI  sequences  in  the 
near  future. 
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results  using  the  ProFit-based  quantitation  of  the  previously 
acquired  GE  1.5T  and  Siemens  3T  and  1.5T  MRS  data  using  2D 
localized-correlation  spectroscopy  (L-COSY)  spectral  acquisition  in 
selected  brain  pathologies  demonstrate  the  increased  power  of 
estimating  more  brain  metabolites  such  as  glutathione  (GSH),  PCh, 
PE,  and  GPC  (109,110).  Shown  in  Figure  11  are  plots  showing  the 
fitting  of  a  2D  JPRESS  prostate  spectrum  (a)  using  a  Siemens  3T 
MRI  scanner,  the  same  as  shown  in  Figure  5,  its  fit  (b)  and  the  fit 
residue  (c)  as  determined  with  the  ProFit  algorithm.  The  estimated 
metabolite  ratios  (/Cr)  using  the  ProFit  quantitation  were  as  follows: 
Cit,  7.36;  PCh,  0.09;  GPC,  0.48;  Spm,  8.55;  ml,  3.79;  si,  0.12;  Glu,  1.52; 
Gin,  0.03;  Tau,  1.47.  These  values  were  in  agreement  with  the  work 
published  by  Lange  et  al.  (69).  However,  the  applicability  of  ProFit 
in  a  clinical  setting  is  yet  to  be  demonstrated. 

FUTURE  DIRECTIONS:  CLINICAL 
APPLICATIONS 

As  discussed  above,  localized  2D  MRS  has  left  infancy  and  moving 
to  maturity.  Twenty  minutes  of  total  duration  were  required  for 
recording  the  decade  old  single-voxel-based  2D  L-COSY  and  2D 
JPRESS  spectra.  Recent  fully  sampled  multivoxel-based  4D  EP-JRESI 
and  EP-COSI  sequences  enable  recording  2D  COSY  and  J-resolved 
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